A topological look at the quantum spin Hall state 
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We propose a topological understanding of the quantum spin Hall state without considering any 
symmetries, and it follows from the gauge invariance that either the energy gap or the spin spectrum 
gap needs to close on the system edges, the former scenario generally resulting in counterpropagating 
gapless edge states. Based upon the Kane-Mele model with a uniform exchange field and a sublattice 
staggered confining potential near the sample boundaries, we demonstrate the existence of such 
gapless edge states and their robust properties in the presence of impurities. These gapless edge 
states are protected by the band topology alone, rather than any symmetries. 
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Since the remarkable discovery of the quantum Hall 
effect (QHE) [l|, the study of edge state physics has 
attracted much attention on both theoretical and ex- 
perimental sides. Recently, a new class of topological 
states of matter has emerged, called the quantum spin 
Hall (QSH) states [3, Q. A QSH state of matter has a 
bulk energy gap separating the valence and conduction 
bands and a pair of spin-filtered gapless edge states on 
the boundary. The QSH effect was first predicted in two- 
dimensional (2D) models @, 13], and was experimentally 
confirmed soon after in mercury telluride quantum wells 
[H • The QSH systems are 2D topological insulators 0, @| 
protected by the time-reversal symmetry (TRS), whose 
edge states are robust against perturbations such as non- 
magnetic disorder. 

A simple model of the QSH systems is the Kane-Mclc 
model |2j, defined on a honeycomb lattice, first intro- 
duced for graphene with spin-orbit couplings (SOCs). It 
was suggested [2j that the intrinsic SOC in graphene 
would open a band gap in the bulk and also establish 
spin-filtered edge states that traverse the band gap, giv- 
ing rise to the QSH effect. Even though the intrinsic SOC 
strength in pure graphene is too small to produce an ob- 
servable effect under realistic conditions Q, the Kane- 
Mele model captures the essential physics of the QSH 
state with nontrivial band topology [3 H . In the pres- 
ence of the Rashba SOC and an exchange field, the Kane- 
Mele model enters a TRS-broken QSH phas e II Oil charac- 
terized by nonzero spin Chern numbers ljj, [l2j . Prodan 
proved [12j that the spin-Chern numbers are topological 
invariants, as long as the energy gap and the spectrum 
gap of the projected spin operator Pa z P stay open in the 
bulk, where P is the projection operator onto the sub- 
space of the occupied bands and a z the Pauli matrix for 



the electron spin. Unlike the Z2 invariant 13[ , the robust 
properties of the spin-Chern numbers remain unchanged 
when the TRS is broken [13, El- 

The existence of counterpropagating edge states with 
opposite spin polarizations is an important characteristic 
of the QSH state. It is believed that the edge states can 
be gapless only if the TRS Q or other symmetries, such 
as the inversion symmetry [l4j or charge-conjugation 



TRS [1511 , are present. When the TRS is broken, it was 
found |10j that a small gap appears in the spectrum of 
the edge states, which was obtained for a ribbon geome- 
try under ideal boundaries, i.e., boundaries created by an 
infinite hard-wall confining potential. However, since the 
edge states are localized around the sample boundaries, 
they can be sensitive to the variation of on-site potentials 
near the boundaries flij . 

In this Letter, in order to reveal the general character- 
istics of the edge states and their connection to the bulk 
topological invariant in a QSH system, we present a topo- 
logical argument similar to the Laughlin's Gedanken ex- 
periment without considering any symmetries. We show 
that, as required by the nontrivial band topology and 
gauge invariance, either the energy gap or the spin spec- 
trum gap (the gap in the spectrum of Po~ z P) needs to 
close on the edges of a QSH system. These two scenarios 
will lead to gapless or gapped edge modes, respectively. 
In particular, it is demonstrated that gapless edge states 
can appear in a TRS-broken QSH system by tuning the 
confining potential at the boundaries. They are associ- 
ated with the bulk topological invariant, and are robust 
against relatively smooth impurity scattering potential. 
Our result offers an interesting example for counterpropa- 
gating gapless edge states that are not protected by sym- 
metries, which sheds light on the underlying mechanism 
of the QSH effect in a broad sense. 

Let us first look back on a looped ribbon of the QHE 
system, with a magnetic flux (f> (in units of flux quantum 
hc/e) threading the ring adiabatically 17 - 19| . The Fermi 
energy Ep is assumed to lie in an energy gap. In the 
spirit of the Laughlin's argument, increasing <j> from 
to 1 effectively pumps one occupied state from one edge 
to the other, giving rise to the transfer of one charge 
between the edges, essentially because there is a nonzero 
Chern number (Hall conductivity) in the bulk. On the 
other hand, the system Hamiltonian is gauge invariant 
under integer flux changes, i.e., if cj> is increased from 
to 1, the system will reproduce the same eigenstates as 
at (j) — 0. To assure this gauge invariance, there must 
be gapless edge modes on the edges (when the perimeter 
of the ring is large), so that the spectral flow can form 
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a closed loop, as illustrated in Fig. la, along which the 
electron states can continuously move with changing <f>. 
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FIG. 1: (a) A schematic of the flow of electron states in a 
looped ribbon of the QHE system, with adiabatically increas- 
ing the magnetic flux <j> that threads the ring. The bulk elec- 
tron states below Fermi energy Ef drift from left to right, 
gapless edge modes ascend through Ef on the right edge, 
states above Ef drift from right to left, and gapless edge 
modes descend through Ef on the left edge, forming a closed 
loop, (b) In the first scenario for the QSH system, gapless 
edge modes appear on the edges of the ring, so that electron 
states in each spin sector behave like in the QHE system, but 
those in two spin sectors move in opposite directions, (c) In 
the second scenario for the QSH system, the edge states are 
gapped, whereas the spin spectrum gap closes on the edges. 
In the bulk, electron states in the two spin sectors drift in 
the same way as in (b), but on the edges they join together 
within the valence (conduction) band. 

We now propose a topological understanding of the 
QSH system in terms of the same looped ribbon geom- 
etry. The occupied valence band can be decomposed 
into two spin sectors by using the projected spin oper- 
ator (The unoccupied conduction band can be 
divided similarly.) The two spin sectors are separated by 
a nonzero spin spectrum gap in the bulk. They carry op- 
posite spin Chern numbers, so that increasing (f> pumps a 
state of the spin up sector in the occupied band from one 
edge to the other, and pumps another state of the spin 
down sector in the opposite direction. In order for the 
system to recover the initial eigenstates as <f> changes from 
to 1, the spectral flow needs to form closed loops, simi- 
larly to the QHE system. However, for the QSH system, 
if not enforced by any symmetry, two different scenarios 
can occur on the edges. One is that gapless edge modes 
appear on the edges, so that states can move between the 
conduction and valence bands with changing cj> to form 
closed loops in the spin-up and spin-down sectors sepa- 
rately, as shown in Fig. lb. In this case, the states in the 
two loops cannot evolve into each other due to the non- 
vanishing spin spectrum gap both in the bulk and on the 
edges. The other scenario is that a closed loop of spec- 
tral flow is formed between the two spin sectors within 



the valence (or conduction) band, as shown in Fig. lc. 
In this case, the spin spectrum gap must vanish on the 
edges, but the energy gap may remain open in the edge 
state spectrum. 




FIG. 2: (a) A schematic of an armchair honeycomb lattice 
ribbon with atom sites in two sublattices being labeled by 
solid dots and hollow dots, where a is the distance between 
nearest-neighbor sites, (b) Profiles of \Vi\ as functions of Xi, 
for an abrupt confining potential (dashed line) and a relatively 
smooth confining potential (solid line). 

The above topological discussion on the QSH system is 
very general, independent of any symmetries. To demon- 
strate the two scenarios in Figs, lb and lc, in what follows 
we take the Kane-Mele model Q for a honeycomb lattice 
ribbon as an example, by taking into account different 
confining potentials near the edges of the ribbon. It was 
shown that in a suitable parameter range, the Kane-Mele 
system is in the QSH phase protected by the TRS, and it 
can become a TRS- broken QSH phase [id ], when a spin- 
splitting exchange field is applied. Consider an armchair 
ribbon along the y direction, as shown in Fig. 2a, in- 
cluding N x dimer lines across the ribbon. (Results for a 
zigzag ribbon are similar.) The boundaries are at x± = 
and Xff x = ^{N x — 1), where the distance between 
nearest-neighbor sites is chosen as the unit of length. The 
Hamiltonian can be written as H = Hkm + He + He 
with 



(ij) v «u» 

+ iV R ^2c\e z ■ (er x d^Cj, 

(ij) 



(1) 



as the Hamiltonian of the Kane-Mele model. Here, the 
first term is the nearest-neighbor hopping term with c\ = 
(c^,c^) as the electron creation operator on site i and 
the angular bracket in standing for nearest- neighbor 
sites. The second term is the intrinsic SOC with coupling 
strength Vso, where a are the Pauli matrices, i and j are 
two next nearest neighbor sites, k is their unique common 
nearest neighbor, and vector dik points from k to i. The 
third term is the Rashba SOC with coupling strength 
Vr. He = gJ2i c l a z c i stands for a uniform exchange 
field of strength g. He represents a sublattice staggered 
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confining potential, which is given by He — Si^i c l c i 
with 



Vi = ±V [e—r+e 



(2) 



where ± is taken to be positive for sites on sublattice A 
(solid dots) and negative on sublattice B (hollow dots), 
as shown in the Fig. [2^,. In Eq. (2), Vi is strongly x depen- 
dent across the ribbon, equal to ±Vb at the edges (x-i = 
and Xi = ZNx)- It decays exponentially away from the 
edges, with a characteristic length £, as shown in Fig. 
[2b- When the ribbon width is much greater than £, Vi 
essentially vanishes in the middle region of the ribbon. 
Here, we note that in the case of a uniform staggered 
potential Vi — ±Vo {\Vi\ being independent of Xi), it was 
shown [l(| that with increasing |T^|, there is a transition 
from the TRS-brokcn QSH phase to an ordinary insu- 
lator state, where the middle band gap closes and then 
reopens. Therefore, for the confining potential Vi given 
by Eq. (2) with large Vb, the ribbon in Fig. [2] can be 
regarded as a TRS-broken QSH ribbon sandwiched in 
between two trivial band insulators. 
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FIG. 3: (color online) The energy spectrum (a, b) and the 
spectrum of the projected spin Pa z P (c, d) of an armchair 
ribbon for £ = 0.1 (a, c) and £ = 4 (b, d), in which L (R) 
stands for states on the left (right) edge, and f (4.) for the up 
(down) spin polarization. The horizontal arrows in (a) point 
to the small energy gaps in the energy spectrum. At £ = 0.1, 
while the energy spectrum (a) is gapped, the spin spectrum 
(c) is gapless. At £ = 4, the energy spectrum (b) is gapless, 
but the spin spectrum (d) is gapped. 

In order to assure the system in the TRS-broken QSH 
state, we set the parameters V so = O.lt, Vr — O.lt, 
and g — O.lt. The length of the armchair ribbon N y is 
taken to be infinite. The energy spectrum of the ribbon, 
together with the corresponding eigenfunctions (p m (k y ), 
can be numerically obtained by diagonalizing the Hamil- 
tonian for each momentum k y in the y direction. The 
calculated energy spectrum of the armchair ribbon with 
width N x = 240, for the confining potential with Vq = 12t 



fixed and two different decay lengths, is plotted in Fig. 
[3^ and Fig. [3J5. One can see easily that the edge states 
appear as thin lines in the middle bulk band gap of the 
energy spectrum. In Figs. [3^, and [3b, the spin polariza- 
tion of the edge states is labeled with | and I, indicating 
that two spin-filtered channels on each edge flow along 
opposite directions. For the nearly hard-wall confining 
potential of £ = 0.1, the sublattice potential V, is nonzero 
only on the outermost armchair lines, similar to the as- 
sumption in Ref. (l6| . In this case, two small energy 
gaps are observed in the edge state spectrum shown in 
Fig. [3^, in agreement with the previous observation [loj . 
as a consequence of the broken TRS. With increasing the 
decay length £ of the confining potential, the energy gaps 
of the edge states become smaller and smaller. As £ is 
large enough, e.g., £ = 4 in Fig. [5Jd, interestingly, the 
edge states become gapless. 

We now calculate the fcj,-dependent spectrum of pro- 
jected spin operator P<j z P, whose matrix elements are 
given by (if m (k y )\ a z \tp n (k y )} with m and n running over 
all the occupied states. By diagonalizing this matrix, the 
spectrum of the projected spin Pa z P can be obtained. 
For the Kane-Mele model Eq. (1), if Vr — 0, a z com- 
mutes with the Hamiltonian H. Therefore, a z is a good 
quantum number. It follows that the spectrum of Pa z P 
consists of just two values ±1, which are highly degen- 
erate. When the Rashba term is turned on, a z and H 
no longer commute, and the degeneracy is lifted. In this 
case, the spectra of Pa z P between +1 and —1 spread 
towards the origin, but a gap remains for a bulk sample 
if the amplitude of the Rashba term dose not exceed a 
threshold 
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For the ribbon geometry, the situation is more compli- 
cated due to the existence of the edges, and numerical 
calculations are performed to obtain the spin spectrum. 
The calculated spectrum of Pa z P for the same param- 
eters as those in Figs. 3a and 3b is shown in Figs. 3c 
and 3d, which exhibits very interesting behavior. For 
the hard- wall confining potential with £ = 0.1, while the 
energy spectrum of edge states are slightly gapped, with 
increasing k y the spectrum of Pa z P continuously change 
between +1 to —1 without showing any gap, correspond- 
ing to the second scenario shown in Fig. lc. On the 
other hand, for a relatively smooth confining potential 
with £ = 4.0, the energy spectrum of the edge states are 
gapless, but the spectrum of Pa z P displays a big gap, 
and the sudden changes happen at the cross points in 
the energy spectrum of the edge states, corresponding to 
the first scenario shown in Fig. lb. From Fig. 3, it follows 
that as long as the system is in the QSH state, a gapless 
characteristic always appears either in the energy spec- 
trum of edge states or in the spectrum of Pa z P, leading 
to the two types of closed loops for the continuous flow 
of the electron states illustrated in Figs, (lb, lc). 

The result shown in Figs. 3b and 3d, for the relatively 
smooth confining potential, is of particular interest. It 
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indicates that gapless edge states can exist in the TRS- 
broken QSH system, accompanied with a gapped spec- 
trum of Pa z P. Such an interesting behavior can be fur- 
ther understood by the following argument. As long as 
the bulk energy gap does not close, the projected spin op- 
erator Pa z P is exponentially localized in real space with 
a characteristic length about A ~ Tlvf/^-e, where vp is 
the Fermi velocity and Ae is the magnitude of the energy 
gap. [HI For the parameter set used in Fig. 3, A is esti- 
mated to be between 1 to 2 lattice constants. When the 
confining potential Vi is varying relatively slowly in space, 
i.e., £ 3> A, one can find that Pcr z P roughly commutes 
with the confining potential. In this case, the confining 
potential is of no influence on the spectrum of Pa z P. 
Since the spin spectrum has a gap in the bulk [l(3], this 
gap remains to open on the smooth edges, as seen from 
Fig. 3d. As a result, the energy gap has to close due 
to the topological requirement, resulting in gapless edge 
modes, as observed in Fig. 3b. 
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FIG. 4: (color online) Eigenenergies of the edge states as a 
function of the magnetic flux <j> threading the looped geometry 
with size 120 x 60 for (a) l = 1.0, U = t, (b) l = 3.0, U = t, 
and (c) Iq = 3.0, Uq = 3t, where the impurity concentration 
is fixed at 1%. The other parameters are taken to be Vso = 
Vr = g = O.lt, Vo — 12t, and £ = 4. Arrows in (a) indicate 
some of the relatively large energy gaps. 

Finally, we wish to discuss the robustness of the gap- 
less edge states found in the present TRS-broken QSH 
system. We consider a N x x N y sample forming a looped 
geometry as that shown in Fig. 1. Ni nonmagnetic im- 
purities are assumed to be randomly distributed in the 
sample at positions R Q with a — 1, ■■■Ni. An extra 
term Hi = Wic\ci is added to the total Hamiltonian 
H to describe the effect of the impurity scattering, where 
Wi = J2 a U{vi — R Q ) with as the position of the i-th 
atom site. The impurity scattering potential is taken to 
be U{ri — R a ) = (Uq/Iq) exp(— |rj — Rq,|/Zo) with l Q as the 
correlation length and Uq the strength of the scattering 
potential. By inclusion of 1/Zq in the prefactor, the area 
integral of the impurity potential is set to be indepen- 



dent of Iq . Figure 4 shows the evolution of the calculated 
eigenenergies of a 120 x 60 system in the band gap upon 
adiabatic insertion of a magnetic flux <p into the ring, for 
three different impurity scattering potentials. The num- 
ber concentration of the impurities is fixed at 1%. For 
a very short correlation length Iq = 1, for which the im- 
purity potential is nearly uncorrelated from one site to 
another, we see from Fig. 4a that at Uq = t, the energy 
levels of the edge states avoid to cross each other as they 
move close, resulting in small energy gaps in the spec- 
trum, as indicated by the arrows. This level repulsion 
behavior is a signature of the onset of backward scatter- 
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When the characteristic length Iq is increased 
to Iq = 3 with Uq = t fixed, corresponding to a relatively 
smooth impurity scattering potential, all the energy gaps 
vanish, as shown in Fig. 4b. The energy levels move in 
straight lines and continue to cross each other, a clear 



indication of quenching of the backward scattering 11 1 . 
Such a level crossing feature is intact when Uq is increased 
up to 3t for fixed lo = 3, as shown in Fig. 4c. We thus 
conclude that the edge states remain to be robust in the 
presence of relatively smooth impurity scattering poten- 
tial of intermediate strength. This result can be under- 
stood based upon an argument similar to that in the pure 
case. When Iq is greater than the characteristic length A 
of the projected spin operator Pa z P, the impurity scat- 
tering potential nearly commutes with Pa z P, and hence 
does not affect much the spin spectrum gap, so that the 
energy gap needs to close on the edges, which explains 
the level crossing behavior of the edge modes. 

In summary, based upon a general topological argu- 
ment without relying on the TRS or other symmetries, 
we show that in a QSH system either the energy gap or 
the gap in the spectrum of Pa z P needs to close on the 
edges. We find that a TRS-broken QSH system can have 
either gapless or gapped edge states, depending on the 
properties of the confining potential near the boundaries. 
The gapless edge states are protected by the bulk topo- 
logical invariant rather than any symmetries, which can 
remain to be robust in the presence of impurities. 
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